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MiniBooNE (Mini Booster Neutrino Experiment) searches for the — > Ue oscillations with Am^ ~ 1 eV'^/c* 
indicated by the LSND experiment. The LSND evidence, when taken with the solar and atmospheric neutrino 
oscillations, suggests new physics beyond the Standard Model. However, this evidence has not been confirmed 
by other experiments. MiniBooNE has completed its first — > i/e oscillation search using a sample of ~1 GeV 
neutrino events obtained with 5.58 X 10^'^ protons delivered to the Booster Neutrino Beamline. The analysis 
finds no significant excess of Ue events in the analysis region of 475 — 3000 MeV. 



1. Introduction 

The long-standing deficits in observed solar elec- 
tron neutrinos (i/g) [D H El H 13 [51 [7] and atmo- 
spheric muon neutrinos (j/^) [51 [Hi UHllIIl [H] now have 
a firm explanation in terms of the flavor oscillations 
of the three neutrinos in the Standard Model. This 
interpretation has been confirmed with experiments 
using terrestial neutrino sources [T31 HH 1151 jTB] and 
point to oscillations with mass-squared differences of 
Am?2~8 X 10-5 eV^c'' and Amig ~ 3 x IQ-^eV^c* 
for the solar and atmospheric oscillations, respectively. 
However, the indications for Vf^ — > Ve oscillations 
reported by the LSND collaboration with Am^ ~ 
0(1 eV^) |T7] cannot be accommodated within this pic- 
ture. Taken with other constraints, such as the num- 
ber of light neutrinos with standard weak couplings 
obtained from measurements of the Z width [18l [19], 
this result would require a dramatic departure from 
the Standard Model in the form of "sterile" neutri- 
nos without standard weak couplings and/or exotic 
forms of symmetry violation [2(7. The LSND results 
have not been confirmed, though other experiments 
searching for neutrino oscillations with similar values 
of Am^ have not had sufficient sensitivity to rule it 
out completely [211 [22] • 

The Mini Booster Neutrino Experiment (Mini- 
BooNE) searches for the neutrino oscillations indi- 
cated by the LSND result using an 0(1 GeV) i'^ 
beam produced by 1.6-/iS-long pulses of 8 GeV pro- 
tons from the Fermilab Booster synchrotron ^23]. Typ- 
ically, '--^ 4 X 10^^ protons are delivered in such a pulse 
at a rate of ~ 4 Hz. The protons in each spill are 
measured by two toroids before they impinge on a 71- 
cm-long beryllium target, where the p-Be interactions 
produce secondary mesons (dominated by pions, with 
a small contribution from kaons) , which in turn decay 
in a 50-meter-long decay region following the target 
to produce the neutrino beam. The beryllium target 
is embedded within an electromagnet ( "horn" ) pulsed 
synchronously with the beam with a 174 kA current 
to produce a toroidal magnetic field that focuses pos- 
itive (negative) particles, resulting in a z/(j7)-enhanced 
beam. In its first results, MiniBooNE has performed 



a search for i/^ — > i/^ oscillations using 5.58 x 10^° 
protons-on-target collected in neutrino mode [24]. The 
resulting neutrino beam is > 99% pure in and peaks 
at ^ 700 MeV, with a small component of coming 
from three-body decays of muons and kaons. 

The MiniBooNE detector [2S| is situated 541 me- 
ters downstream of the target, near the axis defined 
by the beam. This distance, along with the energy 
of the i^^ beam, matches the distance/energy ratio of 
the LSND experiment ('-^ 0.8m/ MeV), with the result 
that the oscillation probability is similar in the two 
experiments. The detector is a 12.2 meter diameter 
sphere filled with 800 tons of undoped Marcol 7 min- 
eral oil with an index of refraction of ~ 1.47. The 
sphere is divided into two concentric, optically decou- 
pled regions at 575 cm radius, resulting in "main" and 
"veto" regions. The main region is instrumented on its 
outer surface with 1280 inward-facing 8" photomulti- 
plier tubes (PMTs) covering 10% of the surface. The 
veto region is instrumented with 240 8" PMTs. Neu- 
trino interactions within the main region are identified 
via the Cherenkov radiation produced by the charged 
particles emerging from the interaction and detected 
on the PMT array. Cosmic muons entering the main 
region can be identified and rejected by the light pro- 
duced and detected in the veto region. The PMT 
activity in a 19.2 fj,s window around each batch of 
protons delivered from the Booster is recorded, where 
the expected arrival time of the neutrinos occurs 4.6 
to 6.2 ^s after the start of the window. Other ac- 
tivity in the detector, in the form of random triggers 
and calibration data, is recorded to study the detector 
response and non-beam backgrounds. 



2. The Booster Neutrino Beam 

The predicted neutrino flux at the detector 
is obtained from a GEANT4-based Monte Carlo 
simulation [25]. The simulation includes a detailed de- 
scription of the beamline geometry, the spatial and 
kinematic properties of the primary proton beam inci- 
dent of the target, and pion and kaon production mod- 
els tuned to the available data. The primary protons 
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Figure 1: Absolute particle production rates in 8GeV p-Be interactions for tt"^ (left) and K'^ (right). The tt^ measure- 
ments show the double differential cross sections measured by the HARP experiment as a function of tt^ momentum in 
bins of the polar angle of the tt"*" relative to the incoming proton. The K'^ measurements show invariant cross sections 
measured with proton energies of 9.5-24 GeV and extrapolated to 8 GeV using the Feynman scaling hypothesis. The 
solid red(black) lines for tt^{K'^) show the parametrization obtained by fitting the data, while the dashed lines show the 
variation in the parametrized cross sections when the parameters are varied according to their uncertainties. 



and secondary particles are tracked through the beam- 
line geometry, accounting for the magnetic field due to 
the horn, as well as electromagnetic and hadronic in- 
teractions in the material (which can produce further 
particles), until a particle decays to produce a neu- 
trino. In the case where a decay chain has multiple op- 
portunities to produce neutrinos {e.g. tt^ — > /i"*" -I- f^, 
/i+ e'^ + V + Vf.) , all decays producing neutrinos 
are recorded. 

The primary source of the neutrino flux is the decay 
of TT"*" from p-Be interactions. The HARP experiment 
at CERN has measured the momentum and angular 
distribution of tt^ production for 8 GeV protons on 
a thin beryllium target [i27^. Similar measurements at 
6 GeV/c and 12.3GeV/c proton momentum are avail- 
able from the BNL E910 experiment. These measure- 
ments are summarized as absolute double differential 
cross sections in pion momentum and angle (relative 
to the incident proton) and parametrized using the 
Sanford-Wang function |28j. The parameters of the 
function are obtained via a fit to the three datasets. 
The results are shown on the left in Figure [T] where 
the parametrization (in red) is overlayed on the 7r+ 
double differential cross sections from the HARP ex- 
periment. The cyan curves show the variations in the 
function when the underlying parameters are varied 
according to the uncertainties returned by the fit. The 
function is used directly in the simulation to deter- 
mine the multiplicity and kinematics of 7r+ mesons 
emerging from p-Be interactions. A similar procedure 
is used for the 7r~ production based on the HARP and 
E910 data, while an E910 analysis of if^ production 
is used to obtain a Sanford-Wang parametrization of 



production. 

The production in p-Be interactions has been 
measured by a number of experiments. However, no 
measurements at 8 GeV proton energy exist. As a re- 
sult, measurements at energies between 9.5 and 24 
GeV[5S] are used, with the Lorentz-invariant cross 
sections extrapolated to 8 GeV using an assumption 
of scaling in the Feynman x variable |30j. The scaled 
invariant cross sections are parametrized as a function 
of X, with the results shown on the right in Figure [T] 
As with the n'^ Sanford-Wang parametrization, this 
Feynman scaling-based function is used directly in the 
simuation to determine the multiplicity and kinemat- 
ics of the outgoing in p-Bc interactions. 

Cross sections for three categories of hadron- nuclear 
interactions (elastic, quasi-elastic, and inelastic scat- 
tering) have been adjusted to match existing measure- 
ments for proton, neutron and charged pion on beryl- 
hum and aluminum|2Il [Ml ISSl IS]- In some cases, 
theoretical guidance [3S] is used to infer cross sec- 
tions where direct measurements do not exist. In all 
other hadron-nucleus interactions, the GEANT4 de- 
fault cross sections are used. 

The predicted neutrino flux at the MiniBooNE de- 
tector is shown in Figure [2] The dominant part of 
the flux, particularly at < 2 GeV, is due to pion 
decays, while the kaon contribution becomes larger 
for energies > 2.5 GeV. For the flux, muon decay 
is the largest component, while kaon decays become 
dominant at > 1 GeV. A number of systematic un- 
certainties in the flux prediction are considered. The 
largest uncertainty arises from the modeling of the sec- 
ondary particle production, while uncertainties in the 
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Figure 2: The neutrino flux at the MiniBooNE detector predicted by a GEANT4-based Monte Carlo simulation of the 
Booster Neutrino Beamline. 



hadronic cross sections also significant. The z/^ flux 
at high energy is particularly sensitive to the details 
of the modeling of the horn magnetic field as well as 
the hadronic cross sections. While systematic uncer- 
tainties in the predicted neutrino fiux, together with 
neutrino interaction cross sections, lead to large uncer- 
tainties in the predicted rate of both background and 
signal processes, in situ measurements of the event 
rates of a number of neutrino-induced processes (dis- 
cussed in Section [s]) reduce the impact of these uncer- 
tainties. 



3. The MiniBooNE detector 

The primary means of (charged) particle detection 
in MiniBooNE is via the Cherenkov radiation pro- 
duced by such particles as they traverse the mineral 
oil. A full understanding of the detector response to 
such particles also requires understanding the scin- 
tillation processes which constitute another source 
of light production, as well as the various processes 
which the optical photons from both sources undergo 
prior to their detection by the PMTs. A detailed un- 
derstanding of the properties of the charge and time 
response of the PMTs to these photons is also needed. 
The measurements and constraints described in the 
following sections are incorporated into a GEANT3- 



based Monte Carlo simulation |36l of the detector. 



3.1. External Measurements on Marcol 7 

While the Marcol 7 mineral oil used in MiniBooNE 
is remarkably transparent (extinction lengths of over 
30 meters at 400 nm, near the peak of the photo- 
cathode sensitivity of the PMTs), it also exhibits a 
rich array of optical phenomena. In addition to the 
Cherenkov and scintillation processes which gener- 
ate light, these include fluorescence processes with 
different excitation/emission spectra and lifetimes, 
Rayleigh and Raman scattering, and absorption. The 
left plot in Figure [3] summarizes the rate of the various 
processes as a function of wavelength. The cumula- 
tive extinction rate (the sum of all optical processes) 
is shown as the black line. In the near ultraviolet re- 
gion ( < 320 nm) , a number of fluorescence processes 
dominate, leading to a large increase in the extinction 
rate in this region. In the visible region (> 320 nm), 
the dominant processes are Rayleigh scattering and 
absorption. 

The optical processes are identified and studied 
in ex situ studies using small samples of Marcol 7 
with 1 — 10 cm path lengths [3T|. The index of re- 
fraction, which summarizes the dielectric properties 
of the mineral oil relevant for Cherenkov radiation 
and Rayleigh scattering, has been measured using 



fpcp07_122 



4 



Flavor Physics and CP Violation Conference, Bled, 2007 




Figure 3: Left: Rates of various optical processes in Marcol 7 from external measurements as a function of wavelength. 
The solid black line is the overall extinction rates obtained from spectrophotometer measurements through a 10 cm cell. 
The dashed black line is the extrapolated extinction rate based on in situ data. The curves labled "Fluor 1-4" are the 
measured excitation rates for the four identified fluorescence processes. The lavender points represent the measured 
rate of Rayleigh scattering at 442 and 538 nm, with the dashed lavender and gray lines representing the theoretically 
extrapolated rates. Right: The measured angular distribution of Rayleigh scattering at 442 nm in the plane perpendicular 
(red) and parallel (blue) to the incident polarization. The latter allows the extraction of the contributions from isotropic 
and anisotropic fluctuations. 



a sodium lamp (589.3 nm) and the observed disper- 
sion used to parametrize the wavelength dependence. 
The extrapolated wavelength dependence is verified 
with direct measurements at wavelengths between 
400 — 700 nm. The wavelength dependence of the 
overall extinction rate is measured by analyzing the 
transmission rate through Marcol 7 samples using a 
spectrophotometer . 

The angular dependence of scattering in Marcol 7 
is determined by measuring the intensity of light scat- 
tered from 442 and 532 nm lasers as a function of 
angle using a PMT as shown on the right in Figure 
|3] Together with its dependence on the polarization 
of the incident and scattered light, the scattering is 
found to be consistent with Rayleigh scattering from 
both isotropic and anisotropic thermal density fluctu- 
ations. The absolute rate of the scattering was deter- 
mined by calibrating the measured PMT rates using 
suspensions of 50 nm polysterene spheres with known 
number density using the scattering cross section cal- 
culated from Mie theory 

The excitation and emission spectra of the fluo- 
rescence properties are determined using steady-state 
spectrophotometer measurements at 250 — 600 nm 
wavelength on 1 cm samples of Marcol 7. A singular- 
value decomposition analysis identified four major 
components of fluorescence. Time-resolved measure- 
ments, which use a pulsed dye laser emitting ultravio- 
let light at wavelengths between 285 and 310 nm with 



7ps width, decomposed the observed time distribution 
of emitted light into components with different emis- 
sion lifetimes as a function of emission wavelength. 
The result is the emission spectrum for each iden- 
tified lifetime component. By matching these emis- 
sion spectra to those determined from the steady-state 
measurements, a complete model of each fluorescence 
process (excitation spectrum, and spectrum and life- 
time of emission) is obtained. Finally, the steady state 
measurements also identified a Raman scattering pro- 
cess. 

3.2. Photomultiplier Response 

The 1520 8" PMTs in MiniBooNE are of two 
types: 322 model R5912 PMTs and 1198 model R1408 
PMTs, both from Hamamatsu. All of the R5912 
PMTs are located in the main array, while the veto 
array is composed entirely of R1408 PMTs. A de- 
tailed understanding of the PMT response to optical 
photons is needed to represent the response of the 
detector to neutrino events accurately. The time re- 
sponse of the PMT plays a critical role in the tuning 
of the optical properties of Marcol 7 using in situ data 
described in Section 13.31 

The response of the PMTs to optical photons is 
characterized by both external and in situ calibration 
data. The wavelength-dependent efficiency of the pho- 
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Figure 4: Left: Reconstructed photon arrival times for R1408 PMTs (top) and R5912 PMTs (bottom) for light flashed 
from the center laser flask. The black histogram is the distribution from data , while the red represents the Monte 
Carlo simulation with reflections and scattering suppressed to isolate the time structure from the PMT response. A flat 
background level has been subtracted from the data distribution based on the observed rate before the expected arrival 
time of the photons. The green and blue histograms the effects of turning on the scattering and reflections, respectively 
so that the latter includes all known optical processes. Right: Correlations between the reconstructed time and charge 
observed in center laser flask data observed in R5912 PMTs. The four bottom plots show the reconstructed charged 
distribution observed in the four regions of reconstructed photon arrival times indicated in the top plot. 



tocathode has been measured by the manufacturer. 
The time response for a small sample of PMTs, partic- 
ularly the late-pulsing behavior, along with the vari- 
ation of the efficiency with incident angle, has been 
characterized in external measurements using a pulsed 
LED [39 . 

Within the MiniBooNE detector, the response of 
the PMTs is studied using a laser calibration system. 
The system consists of four glass flasks situated within 
the main region (including one at the center) contain- 
ing colloidal silica (Ludox@[30]) to disperse light en- 
tering the flasks via optical fibers coupled to pulsed 
diode lasers at 397 and 438 nm wavelength outside 
the detector. The primary purpose of the system is 
to provide gain and time-offset calibrations via the 
charges measured in low-intensity pulses and the re- 
constructed times from the center flask. The low in- 
tensity minimizes photon pileup so that the pulses 
recorded by the PMTs correspond to single photo- 
electrons. Due to the controlled nature of both the 
geometry and timing of the photon emission, the cen- 
ter laser flask system proved to be a valuable tool for 
analyzing the details of the PMT response. Among 
these features include the tube-to-tube variations in 
the delay of the late pulsing and the correlations be- 
tween the recorded charge and reconstructed time rel- 
ative to the expected arrival time. 

The left plot on Figure |4] shows the arrival time of 



hits recorded in laser events relative to their expected 
arrival time. The red, green and blue histograms are 
the results from the Monte Carlo simulation of these 
events, incorporating known PMT and optical effects. 
In the red, optical effects such as reflections and scat- 
tering are shut off, so that the time structure is due 
to effects in the PMT response, such as resolution 
and pre- and late-pulsing. In the green and blue his- 
tograms, the effects of scattering and reflections (off 
the wall and the PMT faces) are respectively turned 
on, resulting in photons which actually have delayed 
arrival times at the PMTs. The right plot shows the 
charge response of R5912 PMTs in different regions 
of reconstructed arrival times in laser events. The 
observed correlations between the charge and time re- 
sponse are modeled within the Monte Carlo simula- 
tion. 

3.3. Tuning with Electrons 

The results of the external measurements form the 
foundation for a model for transporting optical pho- 
tons through Marcol 7 to their detection by the PMTs 
in the detector Monte Carlo simulation. The pa- 
rameters of this model are further tuned using elec- 
trons from cosmic muons decaying-at-rest (/^-DAR) in 
the detector. These events provide a large sample of 
easily-identifiable electron events with well-known en- 
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ergy spectrum that can be reconstructed (as described 
in Section |5| and simulated within the detector Monte 
Carlo. The parameters are constrained such that the 
electrons from /x-DAR simulated using these parame- 
ters are consistent with those observed in data accord- 
ing to a number of quantities such as the angular and 
time distribution of the light measured by the PMTs. 

The process starts by randomly drawing a set of 
parameter values that are consistent with the uncer- 
tainties in the parameters. The starting uncertainties 
are based in large part on the errors from the exter- 
nal measurements and physical boundaries. A sample 
of electrons from /i-DAR is simulated with the Monte 
Carlo using these parameters and reconstructed as in 
data to determine the position, direction and energy. 
Based on these reconstructed quantities, the distri- 
bution of some target variable (typically energy, pho- 
ton arrival times or the geometric distribution of the 
detected photons) is produced. This target distribu- 
tion is compared with the corresponding distribution 
obtained from the data and the between the two 
calculated. 

In each iteration of the procedure, a number of such 
parameter vectors are drawn, and a sample of elec- 
trons is simulated and reconstructed with each con- 
figuration to obtain the target distributions. Based 
on the probability from the comparison of the tar- 
get distributions with the data, a weight is assigned 
to each vector, resulting in larger weights for parame- 
ter vectors for which the target distributions matched 
the data well, and small weights for those in which the 
match was poor. The set of weighted vectors are used 
to produce a new space of allowed vectors that incor- 
porates the constraint of requiring the target distribu- 
tion to agree with the data, resulting in an effectively 
smaller space. New vectors can be drawn from this re- 
duced space and the process is iterated with additional 
target distributions, higher statistics, or localizing the 
comparison to different classes of events. Examples of 
the latter include a comparison of the energy distri- 
bution in bins of reconstructed radius of the electrons 
from the center of the detector. 

The determination of the rate of fluorescence and 
scintillation processes, both characterized by delayed, 
isotropic light, is particularly difficult with only elec- 
trons, which are relativistic for nearly their entire path 
length through the detector regardless of energy. In 
order to provide an additional handle, neutrino neu- 
tral current elastic (NCEL) scattering events are used. 
These events produce a recoil nucleon in the event 
with the outgoing neutrino leaving no trace in the de- 
tector. Since the nucleon is typically below Cherenkov 
threshold, the NCEL sample allows one to study scin- 
tillation in an environment where Cherenkov-induced 
fluorescence is suppressed. 

The result of the process is a space of optical param- 
eter values in which /^-DAR electrons simulated with 
parameter values drawn from the space will have prop- 



erties which match those obtained in data. The al- 
lowed space is significantly constrained from its initial 
configuration based on external measurements. The 
constrained space is the foundation for the estimate 
of systematic uncertainties due to detector response 
in the analysis. Neutrino samples representing both 
the signal and background processes are simulated us- 
ing a number of parameter vectors drawn from the 
space. The systematic uncertainties in a given distri- 
bution are obtained by considering the covariance in 
the distribution across the systematically varied sam- 
ples. The central values of the parameters are ob- 
tained separately by a more ad hoc process. In cal- 
culating the covariance for a given distribution, the 
covariance is taken about the distribution from the 
central value, rather than the mean of the distribu- 
tions obtained from the systematically varied Monte 
Carlo samples. 



4. Signal and Background Properties 

At 0(1 GeV) energies, the dominant contribution 
of charged current (CC) neutrino interactions is from 
the charged current quasi-elastic (CCQE) scattering 
process: 

lye + n^e- +p (1) 

At MiniBooNE, approximately 40% of all neutrino in- 
teractions in the detector are i'^ CCQE interactions. 
The CCQE process is a particularly attractive signal 
process for Cherenkov detectors in that the proton is 
typically below Cherenkov threshold, so that the light 
produced in the event is dominated by the outgoing 
lepton, resulting in a simple event topology consist- 
ing of a single Cherenkov ring. The signal signature is 
therefore a Cherenkov ring from the electron emerging 
from a i^e CCQE interaction. Furthermore, the quasi 
two-body kinematics of the process allow the determi- 
nation of the incident neutrino energy, up to the initial 
motion of the target nucleon in the nucleus, if the mo- 
mentum of the outgoing lepton is reconstructed. We 
refer to the neutrino energy inferred from the lepton 
energy and direction under the assumption of CCQE 
kinematics as E^^. 

The backgrounds come from misidentified inter- 
actions and the "intrinsic" Ve present in the neutrino 
flux from muon and kaon decays in the beamline. The 
dominant component of misidentified j/^ events arises 
from high energy photons, either from neutral current 
(NC) 7r° production (mainly from neutrino-induced 
A resonances) where the tt" decays to produce two 
photons, or from A — > A^7 decays (also from reso- 
nant neutrino scattering) . These photons induce elec- 
tromagnetic showers much like electrons, resulting in 
Cherenkov rings which mimic the signal Vf, CCQE pro- 
cess. The NC 7r° background can be suppressed to 
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Figure 5: Left: Time distribution of PMT hits observed in a neutrino interaction. The first group of hits at ~ 5.9 /is 
is consistent with a neutrino interaction in the region occuring in the beam arrival window of 4.6-6.2 fis. The second, 
smaller group of hits at ~ 8.4 fis is consistent with an electron from the decay-at-rest of a muon produced in the neutrino 
interaction. Right: Subevent times observed in the beam data. The time on the horizontal axis is such that the start of 
the beam arrival is at ns. The excess due to beam-induced events can be seen when all subevents are considered (black). 
The rising structure within the beam arrival window and the tail after it are due to electrons from the decay-at-rest 
of muons produced in neutrino interactions. Electrons from /i-DAR, both from cosmic and neutrino-induced muons, 
are eliminated by requiring > 200 main hits (red). Cosmic muons are suppressed by requiring < 6 veto hits (green). 
Requiring both > 200 main hits and < 6 veto hits results in a neutrino sample with negligible cosmic background (blue), 
indicated by the absence of subevents occurring outside the beam arrival window. 



the extent that the second photon from the 7r° de- 
cay can be detected, while there are no obvious han- 
dles for suppressing the radiative A background. Pho- 
ton backgrounds also arise from neutrino interactions 
outside of the detector ( "dirt events" ) . The intrinsic 
background is inherently irreducible, as it is identical 
to the signal process in every way apart from their en- 
ergy spectrum. Constraints and cross-checks on both 
sources of background will be considered in Sections 
in and 1 

The NUANCE event generator package [H] is used 
to simulate the neutrino interactions. Based on the 
predicted flux and the elemental composition of the 
mineral oil, NUANCE will assign an interaction chan- 
nel and the final state particle configuration. The 
effects of interactions within the nucleus, such as 
the charge exchange and absorption of pions and the 
A + — > N + N process, are simulated according 
to a model of the carbon nucleus, as are the emission 
of photons from the de-excitation of the nuclear rem- 
nant. The neutrino cross sections are tuned to the 
available data. The parameters of the CCQE process 
are adjusted using the observed distribution in the 
the CCQE sample described in Section [8j 

5. Event Reconstruction 
5.1. Subevent Identification 

The reconstruction of beam trigger events proceeds 
by analyzing the time structure of the PMT hits 



recorded during the 19.2 /is window surrounding the 
expected arrival time of the neutrino beam. Clusters 
of hits in time, referred to as subevents, that exceed a 
threshold 10 PMT hits are identified. Neutrino inter- 
actions will produce a subevent 4.6 — 6.2 /xs after the 
start of the window, when the beam is expected to ar- 
rive. They may produce additional subevents if muons 
are produced and come to rest within the main region 
of the detector, leading to decay electrons some mi- 
croseconds later. The left plot of Figure [s] illustrates 
the time distribution of PMT hits observed in such 
an event. The first group of PMT hits is consistent 
in time with a neutrino interaction from the arriving 
beam, while the second, smaller group of hits is consis- 
tent with an electron from fi-DAH of a muon produced 
in the neutrino interaction. These electrons are at 
lower energy (< 52.8 MeV) than the Ve CCQE signal 
events and can be readily distinguished by their low 
multiphcity of main PMT hits (< 200). As illustrated 
above, these electrons are particularly useful in that 
they can tag the presence primary muons produced in 
CC interactions, as well as 7r+ which undergo the 
7r+ 11'^ e+ decay chain. 

A subevent can be identified as either an incoming 
cosmic muon or an uncontained neutrino interaction 
by the number of veto hits it contains: a single pen- 
etration through the veto by a muon typically leaves 
18 veto hits, while a through- going muon that enters 
and exits the detector will typically leave > 25 veto 
hits. With a simple combination of requirements on 
the number of main and veto hits, subevents due to 
cosmic activity can be eliminated, leaving a pure neu- 
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Figure 6: Left: Efficiency of tiie i/^ selection as a function of reconstructed neutrino energy following the preselection. 
The blue, red and black points show the cumulative efficiency of applying the log(yCe/'C^), log{€.e/ jC-n) and m^^ cuts. 
Right: Expected backgrounds as a function of reconstructed neutrino energy. The red and blue components show the 
intrinsic background from kaon and muon decay, respectively, along with the various contributions from misidentified 
i/fi events stacked to give the total background. The black dashed line shows the total distribution expected in the 



presence of 



oscillations with oscillation parameters corresponding to the best fit values from LSND. 



trino sample as illustrated in the right plot of Figure 

5.2. Track Reconstruction 

Following the identification of subevents, the PMT 
hits associated with a subevent are fit under four hy- 
pothesized track configurations for the the outgoing 
particles in the event. The first pair of fits assume a 
single track, either an electron or a muon, emerging 
from a certain point and time (four parameters) in the 
detector, in a certain direction and with a certain ki- 
netic energy (three parameters) . Based on a vector of 
parameter values and a track hypothesis, a predictive 
model determines the expected number and arrival 
times of photoelectrons at each of the main PMTs. 
Convoluting these predicted values with a model of 
the PMT response, the actual charge and time of the 
hits in the subevent (including PMTs which did not 
register hits) can be compared to the model predic- 
tion to calculate a likelihood. Using the standard 
method of maximum likelihood, the best-fit param- 
eter values are extracted by varying the seven param- 
eters in such a way as to maximize the likelihood using 
MINUITgH]. We note that the single electron track 
model corresponds to the signal i^e CCQE hypothesis, 
while the single muon track model corresponds to the 
the Vfj, CCQE hypothesis. Both fit models assume that 
the outgoing lepton is the dominant source of light in 
the event, with negligible contribution from the recoil 
hadron system. 

A second set of fits involve a twelve parameter 



model with two electron tracks. In this model, two 
tracks emerge from positions displaced from the ver- 
tex, with each track pointing back to the vertex. This 
models the 7r° — > 77 decay, with the displacement 
of each of the tracks from the vertex accounting for 
the distance traversed by each photon before shower- 
ing. The model assumes that the detector response to 
each photon is the same as an electron at the conver- 
sion point with the same time, direction and energy. 
The twelve parameters correspond to the position and 
time of the vertex (four parameters), and the direc- 
tion, energy, and displacement (four parameters) of 
each track. 

A second two-track (2T) fit involves the same model 
with the kinematics of the two tracks confined in such 
a way that the invariant mass is equal to the nom- 
inal TT*^ mass, resulting in a model with eleven free 
parameters. The 2T fits, which use the results of the 
one-track electron fit as a starting point, proceed con- 
ceptually in the same way as the one-track fit in that 
the model produces the predicted time and charge at 
each PMT in the subevent based on the configuration 
of the two tracks, and employ the maximum likelihood 
method to extract the best-fit parameter values. Prac- 
tically, the 2T fits differ in that a number of different 
configurations are tested for each subevent, with the 
most promising configurations used as starting con- 
figurations for the likelihood maximization process in 
MINUIT. The results of the fits are a complete kine- 
matic reconstruction of the tt'^ — > 77 decay, with and 
without a mass constraint. In particular, with the 
momentum of each photon reconstructed, one can in- 
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fer the invariant mass of the 77 system (trivial in the 
case of the fixed-mass 2T fit) and the momentum of 
the TfO. 



6. Event Selection 

The selection of candidate ly,, CCQE events from the 
beam data starts with a set of simple requirements to 
eliminate obvious Vfj, CC events (only one subevent 
identified in the event), incoming cosmic ray muons 
(less than six veto hits in the subevent) and electrons 
from /i-DAR (greater than two hundred main hits in 
the subevent). The average time of main PMT hits 
in the subevent is required to be consistent with the 
expected beam arrival window. 

This selection identifies a sample of neutrino can- 
didates with neglible contamination from cosmic 
backgrounds that then undergo track reconstruction. 
Based on the electron fit, the reconstructed vertex of 
the event is required to be within 500 cm of the center 
of the detector, where the the radius defined by the 
face of the main PMTs is 548.5 cm. The results of 
the muon fit are used to determine the projected end- 
point of the track assuming that it is a muon, which 
is required to be within 488 cm of the center of the 
detector. This eliminates CC events in which the 
muon decays close to or behind the main PMT array, 
resulting in a decay electron that may evade detec- 
tion. The reconstructed energy from the electron fit 
is required to be greater than 140 MeV, while the re- 
constructed neutrino energy (E^^) is required to be 
in the range 475 - 3000 MeV. 

A number of criteria are imposed to reduce 
the background from i^^ interactions. A cut on 
log{Ce/Cf^), the logarithm of the ratio of the likeli- 
hoods returned from the electron and muon fits, re- 
duces the remaining i/^t CC background. This quan- 
tity tests whether the event fit better under the elec- 
tron track hypothesis or the muon-track hypothesis. 
Similarly, a cut on log{Ce/ Ct^), where C-^ is the likeli- 
hood returned by the fixed-mass 2T fit, rejects events 
which fit better under the tt" hypothesis than the elec- 
tron hypothesis. Finally, a cut on events which have 
a large reconstructed mass from the free-mass 2T fit 
are also rejected. In all three cases, the specific cut 
values are a function of the reconstructed energy from 
the electron fit and optimized simultaneously using 
Monte Carlo-simulated samples to maximize the sig- 
nal sensitivity. 

A blind analysis is implemented by sequestering 
from study any data events which pass all the selection 
criteria (including the number of such events) until all 
aspects of the analysis are finalized. The unblinding 
procedure is described in Section [TO] In the process of 
performing cross checks on background rates, two sub- 
samples of the signal events which are dominated by 



particular backgrounds were unblinded as described 
in Section [9l 

Figure [6] shows the signal efficiency for i^e CCQE 
events for the three t'^ cuts as a function of E^^ and 
the expected background distribution in E^^ . The 
denominator for the efficiency plot is all events in the 
E^^ bin satisfying the main and veto hit cuts as well 
as the fiducial requirements. The background plot is 
normalized to show the expected background rate for 
the 5.58 X lO^*' protons-on-target used in this analysis. 
This plot also shows the expected signal distribution 
for the LSND best-fit oscillation parameters. Table |l] 
summarizes the expected signal and background rates 
for events with 475 < E^^ < 1250 MeV, where most 
of a potential oscillation signal is expected to lie. 

Table I Expected background rates and systematic un- 
certainties for the CCQE selection in the energy range 
475 < < 1250 MeV. Tlie expected signal yield for 

a 0.26% — > Ve transmutation is also shown. The sys- 
tematic uncertainties are correlated across different back- 
ground sources in some cases. 



Process 


Number of Events 


Uf, CCQE 


10 ± 2 




7± 2 


Miscellaneous Events 


13 ±5 


NC 11° 


62± 10 


NC A ^ A/'7 


20 ±4 


NC Coherent & Radiative 7 


< 1 


Dirt Events 


17 ±3 


from ^ Decay 


132 ± 10 


Ve from Decay 


71 ± 26 


Ve from Decay 


23 ±7 


from TT Decay 


3± 1 


Total Background 


358 ± 35 


0.26% i/p 


163 ± 21 



7. Boosted Decision Tree Selection 

A second event selection system for signal Ve CCQE 
candidates based on a boosted decision tree (BDT) 
|44] is used to cross check the likelihood ratio-based se- 
lection. This method uses a simpler event reconstruc- 
tion with corresponding electron, muon and two-track 
hypotheses. Based on this fit, 172 variables character- 
izing the space and time distribution of the light are 
formed using quantities like the charge and time likeli- 
hoods in different angular regions, the likelihood ratios 
formed from the three fits, and the parameters from 
the two-track fits. These quantities form the input 
for the BDT, which is trained using a "cascade" tech- 
nique f45^ on simulated samples of background and 
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signal events. For each event, the BDT polls the out- 
put of its decision trees to determine the final output. 



8. Background Measurements 

The rate of a number of background processes can 
be measured or inferred from the neutrino data col- 
lected at MiniBooNE. One such measurement is based 
on the identified i/^ CCQE events which provide a 
measure of the dominant i^^ flux. A second measure- 
ment is the rate of NC tt'^ events. In both cases, the 
direct determination of the rate within the data allows 
one to reduce systematic uncertainties associated with 
the Monte Carlo-based prediction of the neutrino flux 
and the neutrino cross sections. 

8.1 . z/^ CCQE events 

The selection of j/^ CCQE events starts by identi- 
fying events with a single subevent consistent with an 
electron from /i-DAR following the primary subevent. 
The primary subevent is required to occur in the beam 
arrival window and have less than six veto hits and 
greater than two hundred main hits in order to sup- 
press cosmic background. The second subevent is re- 
quired to have less than six veto hits and less than 
two hundred main hits, rejecting cosmic muons and 
ensuring that the energy is consistent with an electron 
from fi-DAR. Following the reconstruction of both 
subevents with the single-track reconstruction, the ge- 
ometry of the two subevents is analyzed to determine 
whether the electron vertex is consistent with the end- 
point of the muon. 

This consistency is determined by comparing the 
muon track length determined from the reconstructed 
muon and electron vertices, and the expected muon 
track length based on the reconstructed muon energy. 
For i^f^ CCQE events, where the muon is the domi- 
nant source of light, this is expected to be linear rela- 
tionship due to the minimum-ionizing behavior of the 
muon. This is demonstrated in the top left figure of 
Figure [7] which shows the reconstructed track length 
versus muon kinetic energy for a sample of Monte 
Carlo-simulated i'^ CCQE events. The bottom fig- 
ure shows the same distribution observed in the data. 
The core of the distribution lies on a "line" with slope 
~ 2MeV/cm, but a cluster of events lies below it, 
corresponding to background events where extra par- 
ticles in the hadronic system disrupt the relationship. 
These particles introduce more energy into the event 
while the reconstructed track lengths, based on vertex 
information, are not affected. This places the events 
below the line (more energy per track length). 

A precise form for this line is determined by pro- 
filing the two dimensional distribution observed in 
data and fitting it to a linear function. By using 



the data directly (rather than Monte Carlo-simulated 
CCQE events), systematic uncertainties associated 
with energy scale are reduced. The function summa- 
rizes the predicted energy /track length relationship 
for CCQE events and is used by comparing the 
distance between the projected muon endpoint (us- 
ing the reconstructed energy, direction and the en- 
ergy/track length relationship) and the reconstructed 
electron vertex. The vector difference between these 
two points is projected along the reconstructed muon 
direction resulting in the longitudinal electron ver- 
tex displacement. This distribution of this quantity 
for two subevent i/^ CC candidates is shown on the 
right of Figure [Tj The background, particularly those 
with pions in the final state, lie towards negative val- 
ues (electron vertex closer than would be expected 
based on the reconstructed muon energy) while the 
t'^ CCQE events are centered at zero. The i/^ CCQE 
candidates are selected by requiring that the longitu- 
dinal electron vertex displacement is less than 50 cm. 

The distribution of the selected i/^ CCQE 

events is used in two ways to determine event rates 
of other processes. The first is the rate and spec- 
trum of signal v^, oscillation events for a given set of 
oscillation parameters (sin^ 29, Am^). For small oscil- 
lation probabilities, the rate of CCQE events can 
be translated into the rate of oscillation events since 
they originate in the same flux of neutrinos and inter- 
act via the same CCQE process. Differences between 
and Ve cross sections are considered in the system- 
atic uncertainties. In order to make this translation, 
one must account for the effects of resolution and effi- 
ciencies. The detector Monte Carlo simulation is used 
to infer the energy spectrum of the Vfj^ events from the 
observed E^^ distribution, which is then used to de- 
termine the rate of oscillation events as a function of 
neutrino energy. 

This determination reduces systematic uncertain- 
ties in the predicted signal yield by using the observed 
rate CCQE rate directly to infer the rate of 
CCQE events. The procedure is also used to correct 
the rate of other i/^-induced processes (apart from NC 
7r° and radiative A events), though the systematic un- 
certainties are typically larger when trying to relate 
the rate of processes which interact with intrinsically 
different neutrino cross sections. 

Second, the dominant source of the intrinsic Ve from 
muon decay (the largest single background source) are 
from decays of the same tt"*" which decay to produce 
the i^^ observed in the z/^ CCQE sample. As a result, 
the rate of muons from 7r+ decays in the beamline 
can be inferred from the observed rate of these events 
in the detector. Furthermore, due to the tight kine- 
matic constraints for 7r+ /i+ + i^^ decays to send a 
neutrinos towards the detector (i.e. only decays with 
forward j/^ are observed) , there is a close relationship 
between the observed energy spectrum and the en- 
ergy spectrum of the parent pions. This allows the 
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Figure 7: Left: Reconstructed muon/electron distance versus visible energy for a sample Monte Carlo-simulated ffj, 
CCQE events (top) and data (bottom) where each sample is required to have two subevents consistent with a neutrino 
interaction and an electron from muon decay-at-rest. Right: the longitudinal displacement of the reconstructed electron 
vertex relative to extrapolated muon endpoint as described in the text. The non-i/^ CCQE backgrounds are shown as 
the solid pink histogram. The dashed pink shows the subset of the background which contains no pions in the final state. 



spectrum of this background to be inferred from the 
neutrino data. Residual systematic uncertainties re- 
sult from the detector response uncertainties and neu- 
trino cross section uncertainties which affect the effi- 
ciency of the i^f^ CCQE selection and the background 
rates. Since the and Ve events which result from 
pions of a given energy have different energy spectra 
(the former are from the two-body decay of the , 
while the latter are from the three-body decay of /i+ 
from the two-body decay of tt"*"), there are additional 
uncertainties associated with the energy dependence 
of the neutrino cross sections. 

8.2. NC 7r° events 

A large fraction of the NC 7r° events are readily 
identifiable by the invariant mass determined by the 
free-mass 2T fit. The events in the tt" peak of the in- 
variant mass distribution are used as a control sample 
in which the rate and kinematic properties of the NC 
7r° events are studied. This rate is used to correct the 
rate and spectrum of the NC 7r° background contri- 
bution to the signal fe CCQE sample. While the ex- 
trapolation between these disjoint sample depends on 
the detector Monte Carlo simulation to predict the ef- 
ficiency of selecting NC tt" events, the rate of non-NC 
tt'^ events within the control sample, and the fraction 
of NC tt'^ events that pass the selection criteria, the 
systematic uncertainties incurred by a rate estimation 
using the predicted neutrino flux and cross sections 
are greatly reduced by measuring the production rate 
directly in the data. 

The event selection for the control sample starts by 
requiring a single subevent with greater than 200 hits 



in the main PMT array and less than 6 veto hits. The 
single subevent requirement is used since NC 7r° events 
are not expected to have muons. The remaining CC 
background is suppressed by requiring log(£e/£^) > 
0.05. Electron-like events are eliminated by requiring 
log(£e/'C7r) < (likelihood ratio favoring tt*^ rather 
than e) and 90 < m^yj < 180 MeV/ (reconstructed 
mass consistent with nominal tt^ mass) . The selection 
results in a NC 7r° sample of very high purity (typi- 
cally > 95%) so that systematic uncertainties due to 
non-NC 7r° events are minimized. 

The selected NC tt° candidate events are divided 
according to their reconstructed momentum into nine 
bins ranging from to 1500 MeV/c, as shown in Figure 
[Sj Non-NC 7r° backgrounds are subtracted according 
to the predicted rates from the Monte Carlo simula- 
tion. The selection efficiency and momentum resolu- 
tion predicted by the Monte Carlo simulation are used 
to infer the absolute rate of NC 7r° events in true mo- 
mentum bins prior to the selection. The NC tt" rates 
in the Monte Carlo simulation are corrected so that 
the predicted rates in the control sample match the 
observed rate. The same momentum-dependent cor- 
rections are applied to the NC tt" events that enter as 
background into the signal i/^ CCQE sample. 

The systematic uncertainties in the 7r° rate analysis 
are dominated by the uncertainties in the detector re- 
sponse and the composition of the NC 7r° events. The 
latter uncertainty is associated with the fact that the 
analysis is performed agnostically to which exclusive 
channel (resonant, coherent, etc.) is producing the 
tt". On the other hand, the details of the final state, 
such as the presence of other pions, the energy of the 
recoil nucleon, etc., affect the efficiency of the selec- 
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Figure 8: Observed va^^ distributions in tiie ■k^ control sample in bins of reconstructed tt" momentum, from lowest 
(top left) to highest (bottom right). The black points show the observed data. The blue histograms show the expected 
contribution from non-vr'^ events, while the red histogram shows the total Monte Carlo prediction with systematic errors 
following the correction procedure described in the text. The dashed grey line shows the Monte Carlo prediction prior 
to the correction. 



tion and the reconstructed kinematics. As a result, 
the extracted yield of NC 7r° events is sensitive to the 
underlying composition of exclusive channels. 

The measured rate of NC tt" production also con- 
strains the rate of the A — > A^7 process, since the 
dominant source of NC ifi production is A produc- 
tion and the two processes are related by the elec- 
tromagnetic and hadronic branching ratios of the A. 
The predicted rate of A — > events is corrected 
according to the observed NC tt" rate. This correc- 
tion is also subject to the uncertainties in the channel 
composition of the observed tt*^ events, most notably 
the uncertainty in the fraction of NC tt'' events that 
originate from A production. 



9. Background Cross Checks 

In order to check the Monte Carlo-based back- 
ground predictions for the signal v^, candidate sam- 
ple, a number of control samples and distributions 
in the data are investigated. The control samples 
are constructed to enrich a particular type of back- 
ground. The distribution of geometric, kinematic and 
background-suppression variables are compared be- 
tween data and the Monte Carlo simulation via a 
test including all sources of error and their correla- 
tions to check for consistency. We describe a few of 
these important cross checks here. 

The behavior of misidentified NC tt" events, the 
dominant source of background, was examined us- 
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Figure 9: Left: The Rtwb (distance of the vertex from the wall in the direction opposite to the track direction) 
distribution for low energy candidate events. Right: Reconstructed neutrino energy (E^^) distribution for events 
with E^'^ > 1600 MeV. The data are shown as black points while the Monte Carlo-predicted background is shown in red, 
where the boxes indicate the uncertainties. Both distributions are absolutely normalized to 5.58 x 10'^" protons-on-target. 



ing sideband regions in the two background suppres- 
sion variables log(£e/'CT) and While the over- 
all rate of these events has been corrected to match 
the observed rate within the control sample of well- 
reconstructed NC Ti'^ events, these sidebands cross 
check the behavior of the background with events 
which are signal-like in one of the variable (signal- 
like in both variables would result in events which are 
in the blinded signal region). The rate and kinematic 
properties of the events, as well as the distribution of 
the background suppression variables, are found to be 
consistent between the Monte Carlo simulation, indi- 
cating that this background is properly modeled. 

A second source of background comes from neu- 
trino interactions outside of the detector which pro- 
duce high energy photons (typically from ir^ decay) 
referred to as "dirt" events. The veto region is ap- 
proximately one radiation length in thickness. This 
results in some of these photons evading detection in 
the veto and showering within the main region, yield- 
ing electron-like Cherenkov rings. They are character- 
ized by their proximity to the edge of the detector and 
their relatively low energy (typically 200-300 MeV). A 
control sample enriched in these events is isolated by 
considering events which otherwise pass the Ve selec- 
tion criteria, but have a reconstructed vertex close to 
the wall and low visible energy. To further isolate the 
contribution from these events, one examines the dis- 
tribution of events in the variable Rtwb, which con- 
siders the distance of the vertex from a 540-cm-radius 
sphere concentric with the detector along the direction 
opposite to the reconstructed track direction. Since 
dirt events are inward-directed and near the edge of 
the main region, they concentrate at low Rtwb- The 
distribution of Rtwb in this control sample is shown 
on the left in Figure [9] The predicted contributions 



in the Monte Carlo simulation for the two major com- 
ponents of events in this sample, namely NC 7r° pro- 
duction within the detector and the dirt events, are 
also shown. The distribution terminates at 40 cm due 
to the requirement that the event vertex lie within 
500 cm of the center of the detector. The consistency 
between the Monte Carlo simulation and the data in- 
dicate that the rate of this background is accurately 
modeled in the simulation. 

Events at high energy {E^^ > 1.6 GeV) have a 
large contribution from intrinsic Ve fi'om kaon decay. 
The expected rate of i^e events from oscillations is also 
small in this energy region, as seen in the right plot of 
Figure [g] .The E^^ spectrum of i^e candidates passing 
the selection criteria with E^^ > 1.6 GeV is shown on 
the right in Figure |9] The yield of events in this con- 
trol sample is consistent with the default background 
prediction, though the systematic uncertainties, pri- 
marily due to the uncertainty in the rate of high mo- 
mentum TT*^ production, are large. The ly^, candidates 
events in this energy range are incorporated in to the 
signal extraction procedure described in Section |10| 
thus constraining the rate of intrinsic background. 



10. Signal Extraction and Unblinding 

The i^e signal is extracted by a fit to the E^^ dis- 
tribution of the selected fe candidates. The observed 
energy distribution between 475 and 3000 MeV is sum- 
marized in a histogram with eight bins, and the equiv- 
alent distribution for the expected background and 
signal for a given oscillation parameter set (sin^ 29, 
Am^) from the Monte Carlo simulation is created. 
The uncertainties in the background prediction and 
the signal are summarized in a covariance matrix. The 
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Figure 10: Left: Observed E^^ distribution for selected Ue CCQE candidate events. The expected backgrounds are shown 
as stacked components. The error bars on the predicted background represent the square root of the diagonal components 
of the covariance matrix of the distribution. The dashed purple hue shows the expected E^^ distribution including the 
contribution of oscillation events from the best fit oscillation parameters. Right: The E^^ distribution of the excess 
of Ve CCQE candidate events, obtained by subtracting the predicted background from the observed distribution. The 
black line shows the E^^ distribution of the best fit oscillation parameters, while the red and blue show the predicted 
excess from two other oscillation parameters from the LSND allowed region. 



best-fit sin^ 29) and Am^ values are extracted by vary- 
ing these parameters and minimizing the between 
the observed data and the predicted background plus 
signal contributions. The signal extraction assumes 
that the energy distribution of the signal takes the 
form of neutrino oscillations driven by one Am^. 

The results of the analysis were unblinded in a 
staged procedure. The primary means of checking 
the behavior of variables and distributions within the 
blinded events were a series of tests on geomet- 
ric, kinematic and background suppression variables 
using the total covariance from both statistical and 
systematic errors. Since the x^ tests must allow for 
the presence of an excess of events from oscillation, 
the signal extraction fit is performed and the expected 
contribution from the fitted signal propagated into the 
distributions in addition to the background. The co- 
variance matrix for the distribution used in the x^ 
tests are also updated. This procedure is performed 
without reporting the results of the signal extraction 
fit. 

In the first step, only the from the comparisons 
(with the exception of E^^) are reported without 
showing the data distributions, nor the expectations 
from the best fit. In the second step, the data and 
best-fit distributions in the variables (with the excep- 
tion of E^'^) are shown with the normalization and 
signal/background composition suppressed. In the fi- 
nal stage, the event yield, along with the E^^ distri- 
bution and the best fit parameters are unblinded, com- 
pleting the unblinding the process. While the E'^^ 
distribution is used to obtain the best fit parameters, 
neither the x^ probability of this fit nor the distribu- 
tion are examined until the final stage. 

In the process of carrying out the first step using an 
energy range of 300 < E^'^ < 3000 MeV, it was found 



that the for the visible energy distribution had a 
small probability of 0.01. No other distributions were 
found to be problematic. Since the test accounts for 
the possibility of a signal, the discrepancy could not 
be due to an excess of events consistent with neutrino 
oscillations with a single Am^. A possible discrep- 
ancy was also observed in the log{Ce/ C^^) /rrijj side- 
band samples described in Section [9] in the form of 
an excess of events over the predicted background at 
low energy. This discrepancy, however, was not large 
enough to result in a poor x^ for the energy distribu- 
tions in these control samples. An examination of the 
unsigned deviation in the visible energy distribution 
of the observed data in the signal sample to the best- 
fit distribution confirmed that the poor x^ was likely 
due to a discrepancy at low energy. 

Based on this information, the impact of increasing 
the E^^ threshold was investigated. It was found that 
restricting the analysis to 475 < E^^ < 3000 MeV 
did not impact the sensitivity of the analysis to neu- 
trino oscillations with a single Am^. After deciding 
to implement this increased threshold in the analysis, 
the unblinding procedure was repeated. No problems 
were found in this second iteration, with all distri- 
butions, including the visible energy distribution, re- 
turning reasonable x^ probabilities. As a result, the 
procedure was taken to completion. 



11. Results 



The E^^ distribution of the signal candidates is 



shown on the left in Figure 10 In the analysis region 
of 475 < EQ'^ < 3000 MeV, 380 events are observed, 
where the expected background is 355 ± 19(stat) ± 
35(sys), corresponding to an excess of 0.55 standard 
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Figure 11: Left: Limits at 90% confidence level on the — > Ve oscillation parameters obtained from the MiniBooNE 
data using the one-dimensional raster scan method. The solid black line shows the limit, while the dashed black line 
shows the projected sensitivity. The solid blue line represents the limit obtained from the BDT analysis of the same 
data. Right: Limits on — > oscillations from MiniBooNE, KARMEN2 and Bugey. The cyan and light blue regions 
on each plot are the parameter space compatible at 90% and 99% confidence level with the LSND oscillation evidence, 
respectively. 



deviations over background. The fit yields a value of 
{sm^2e),Am^) = (1.1 x lO'^, 4.1 eV^C*). The back- 
ground and signal distribution corresponding to these 
parameters is shown as a dashed purple line. The dif- 
ference in between the null hypothesis and the best 
fit is 0.94, while the corresponding difference for 
the LSND best fit parameters is 13.7. This indicates 
that the former is an adequate fit to the data, while 
the latter is highly disfavored. In summary, both the 
overall yield of events and the E^^ distribution are 
consistent with the expected background. 

A significant excess of events is observed at 300 < 
E^^ < 475 MeV, below the analysis threshold, as ex- 
pected from the poor x^s in the initial consistency 
tests described in Section 10 The right plot in Fig- 



ure [To] shows the excess of data events over back- 
ground. The black histogram shows the expected ex- 
cess for the best fit parameters, while the red and blue 
histogram show the excess expected for two differ- 
ent oscillation parameters consistent with the LSND 
evidence, (sin^ 26*), Am^) = (0.004, 1.0 eV^c'') and 
(0.2, 0.1 eV^/c"*), respectively. While lower values of 
Am^ produce oscillation signatures which concentrate 
at lower energies, it is not possible to accommodate 
the observed excess with neutrino oscillations driven 
by a single Am?. This also follows from the fact that 
the fit did not find an adequate solution when the low 
energy region was included. The excess at low energy 
is currently under investigation. 

The 90% confidence level limits on the i/n i^e os- 



cillation parameters are shown on the left in Figure 
11 where parameters to the right of the solid black 
line are excluded. The limits are obtained in a one- 
dimensional raster scan whereby a limit on the max- 
imum allowed sin^ 29) is determined at each Am^ to 
obtain the curve. For comparison, the dashed line in- 
dicates the projected sensitivity of the analysis. The 
obtained limit is somewhat worse than the projected 
sensitivity due to the small excess of events observed 
in the data. The limits from the BDT analysis, which 
was unblinded simultaneously and did not observe a 
significant excess of events, are also shown in Figure 
[TT| In addition to a different event selection scheme, 
the BDT analysis also used a different signal extrac- 
tion procedure in which the v^^ and CCQE samples 
were fit simultaneously. 



12. Conclusions 

MiniBooNE has searched for i/^j, v,, oscillations 
in a sample of 5.58 x 10^° protons-on-target deliev- 
ered to the Booster Neutrino Beam in neutrino mode. 
The primary backgrounds to the analysis are con- 
strained with in situ measurements and cross checks. 
The analysis yields no evidence for neutrino oscilla- 
tions: the observed yield and energy spectrum of the 
selected v,, candidates are consistent with the back- 
ground and incompatible with the oscillations indi- 
cated by the LSND experiment. An excess of events at 
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energies below the analysis E^^ threshold of 475 MeV 
remains under investigation. The oscillation parame- 
ters excluded by the analysis are shown on the right 
in Figure [TT] together with the allowed regions from 
the LSND analysis and limits from the KARMEN and 
Bugey experiments. 
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